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One-loop renormalization factors and mixing coefficients of bilinear 
quark operators for improved gluon and quark actions * 
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We calculate one-loop renormalization factors and mixing coefficients of bilinear quark operators for a class of 
gluon actions with six-link loops and 0(a)-improved quark action. The calculation is carried out by evaluating 
on-shell Green's functions of quarks and gluons in the standard perturbation theory. We find a general trend 
that finite parts of one-loop coefficients are reduced approximately by a factor two for the renormalization-group 
improved gluon actions compared with the case of the standard plaquette gluon action. 
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1. Introduction 

In numerical studies of field theories through 
computer simulations finite lattice spacing effects 
poses an obstacle in extracting physical predic- 
tions for observable quantities. In an effort to 
reduce this problem, a recent full QCD simula- 
tion by the CP-PACS CollaborationjjJ employs a 
combination of the 0(a)-improved clover quark 
actional and a renormalization-group- improved 
gluon action with 6- link loops ||, based on pre- 
vious results of a comparative study of various 
action combinations p|. With this development 
it has become necessary to calculate renormaliza- 
tion factors and mixing coefficients of quark oper- 
ators for improved gluon and quark actions. We 
have carried out such a calculation to one-loop 
order of perturbation theory for bilinear quark 
operators |^B| . In this article we summarize the 
results. 

Our calculations are performed within stan- 
dard perturbation theory for on-shell Green's 
functions of quarks and gluons. For the plaque- 
tte gluon action, perturbative calculation of the 
mixing coefficients of bilinear quark operators has 
been previously carried out using the Schrodinger 
functional technique @||. A recent work with the 
same technique ||] extends the results to the same 
set of improved gluon actions as examined by us. 
We find a good agreement between results by us 



and those reported in these references. 

2. Choice of actions 

We adopt the clover action for quark, which is 
given by adding the clover term 



to the ordinary Wilson fermion action. The gluon 
action we consider is defined by 
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Tr {collpi + ciU r tg + c 2 U ch r + c 3 U p ig} , 
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where the first term represents the standard pla- 
quette term, and the remaining terms are six- 
link loops formed by a I x 2 rectangle, a bent 
1x2 rectangle (chair) and a 3-dimensional par- 
allelogram. The coefficients cq, •••,03 satisfy 
co + 8ci + I6C2 + 8C3 = 1. At the one- loop 
level, the choice of the gluon action is specified 
by the pair of numbers C\ and c 23 = c 2 + C3. 
We adopt the following five cases in our calcu- 
lation: (i) the standard plaquette action c\ = 
0, C23 = 0, (ii) the tree-level improved action 
in the Symanzik approach c\ = — I/I2,C23 = 
0(nj, and (iii) three choices suggested by an ap- 
proximate renormalization-group analysis, c\ — 
-0.331, c 23 = and ci = -0.27, c 23 = -0.04 by 
Iwasakig, and c x = -0.252, c 23 = -0.17 by Wil- 
son ra. 
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Table 1 _ 

Finite part zp of renormalization factor for bilinear operator of form (1 + ma)ipTip (i.e., z = in (|l|)) 
with m the subtracted quark mass incorporating one-loop self-energy correction. Coefficients of the term 
Cg W (n — 0, 1,2) in an expansion z-p = z r °^ + cgwZjp + c sw z r ' are gi yen m t ne column labeled (n). 
Errors are at most in the last digit given. 



gauge action 

Cl C 2 3 


V 

(0) (1) (2) 


A 

(0) (1) (2) 


S 

(0) (1) (2) 


P 

(0) (1) (2) 



-1/12 
-0.331 
-0.27 -0.04 
-0.252 -0.17 


-20.618 4.745 0.543 
-16.603 4.228 0.464 
-11.099 3.326 0.336 
-11.540 3.418 0.353 
-10.525 3.248 0.338 


-15.797 -0.248 2.251 
-12.540 -0.198 2.021 
-8.192 -0.125 1.610 
-8.523 -0.131 1.657 
-7.707 -0.117 1.587 


-12.953 -7.738 1.380 
-9.607 -6.836 1.367 
-4.858 -5.301 1.266 
-5.260 -5.454 1.292 
-4.366 -5.166 1.287 


-22.596 2.249 -2.036 
-17.734 2.015 -1.745 
-10.673 1.601 -1.281 
-11.292 1.644 -1.316 
-10.001 1.565 -1.212 



3. Bilinear operators and improvement 

We consider bilinear quark operators of form, 

O r = ipTip, r = 1,75,7^,7^75,^. (1) 

At the tree-level, there exists a one-parameter 
family of 0(a)-improved operators given by[p"2"| 

0o = [1 +ar (1 - z)to ]#Y 

+^r®v - z 2 vtv, (2) 

where and V are 0(a) and 0(a 2 ) vertices, 
and mo is the bare quark mass. 

At the one-loop level, on-shcll Green's func- 
tions of these operators do not have terms of 
0(g 2 a\oga) when we set the clover term coef- 
ficient csw = 1 @- The 0(g 2 a) terms still re- 
main, however. In order to extract the renormal- 
ized operators without 0(g 2 a) errors we need to 
calculate mixing coefficients Br and Cr as well as 
the renormalization factor Zr defined by 

T Q = Z^Ol - g 2 C F am R B T O r R 

-g 2 C F aC r dl, (3) 

where Cf denotes the second-order Casimir 
eigenvalue for the quark field, and the last two 
terms are needed to remove 0(g 2 a) errors from 
on-shell matrix elements, with a dimension 4 
operator with derivative, e.g., (D A — d^O p and 
O v = d^O T . 

4. Strategy of one-loop calculation 

We apply standard perturbation theory to the 
Green's functions of the operator Oq and two on- 



shell external quarks to calculate renormalization 
factors and mixing coefficients. In order to iso- 
late terms of 0(g 2 a), one-loop amplitudes are ex- 
panded in powers of a, or equivalently in the ex- 
ternal momenta p and quark mass m. Infrared 
divergences that generally appear in on-shell am- 
plitudes are regularized by a gluon mass A. 

For this procedure to work, infrared diver- 
gences that contribute to the mixing coefficients 
Br and Cr have to cancel out among diagrams. 
We have explicitly checked that this is in fact the 
case; infrared divergences appear only in the ordi- 
nary renormalization factor Zr, which, however, 
cancel against those in the renormalization fac- 
tor for the continuum operator. As a result, the 
renormalization factor in the MS scheme takes 
the form, 

where /i2(T) is an integer given by h,2(T) = 
4(A), A(V), 16(P), 16(5), 0(T). 

We note that the bare quark mass mo in (Q) 
may be replaced by a subtracted mass m = 
Too — g 2 CpJ^o/a incorporating the one-loop self- 
energy correction by a redefinition of zp, -Br & n d 
Cr . We adopt this definition for numerical results 
presented below. 

4.1. Results 

The calculational procedure described above 
yields z-p, Br and Cr as sums of one-loop inte- 
gral constants. While straightforward in princi- 
ple, the actual algebra to evaluate the integrands 
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Table 2 



Mixing coefficients for axial vector, vector, pseudo scalar and scalar density at z = 0. 



gauge action 


C A 


C v 


B A 


By 


B P 


B s 


Cl 


C23 




















-0.005680(2) 


-0.01226(3) 


0.1141(1) 


0.1150(2) 


0.1148(1) 


0.1444(2) 


-1/12 





-0.00451(1) 


-0.01030(4) 


0.0881(1) 


0.0886(2) 


0.0890(1) 


0.1144(2) 


-0.331 





-0.00285(1) 


-0.00729(4) 


0.0547(1) 


0.0550(2) 


0.0561(1) 


0.0747(2) 


-0.27 


-0.04 


-0.00302(1) 


-0.00757(4) 


0.0572(1) 


0.0575(2) 


0.0586(1) 


0.0777(2) 


-0.252 


-0.17 


-0.00281(1) 


-0.00705(4) 


0.0512(1) 


0.0514(2) 


0.0527(5) 


0.0706(2) 



is extremely tedious. We employ Mathematica 
to carry out this task, and also to generate the 
FORTRAN code from the results. Infrared di- 
vergences that can appear in individual integrals 
are subtracted by working out the leading term 
of integrands for small loop momenta. Integrals 
are then evaluated by the Monte Carlo routine 
VEGAS in double precision. We generally em- 
ploy 20 sets of 10 5 points for integration, except 
for Ca for the plaquette action for which 20 sets 
of 10 6 points are used. Errors are estimated from 
variation of integrated values over the sets. 

We present our results for zr, Cr and Br for the 
case of z = in TablesJD and |^. Results for z =/= 
can be found in Refs. . The Wilson parameter 
is taken to be r — 1. For zr, the coefficients of 



(2) 



t-l. • (0) , (1) 

the expansion zp = z T + cswz T -r ^sw^r 
given (Table |l|), while the constants Cr and Br 
are evaluated for csw = 1 (Table ||) . 

We observe in the results that the renormaliza- 
tion factors and mixing coefficients are reduced by 
about a factor two for renormalization-group im- 
proved gluon actions compared to those for the 
plaquette action. More generally we find that 
zr's monotonically decrease as c\ and C23 become 
large and negative]^]. 

Numerical calculations have been carried out 
at Center for Computational Physics, Univer- 
sity of Tsukuba and at Research Institute for 
Fundamental Physics, Kyoto University. This 
work is supported in part by the Grants-in- 
Aid of the Ministry of Education (Nos. 2373, 
09304029,10640246). KN and YT are JSPS Re- 
search Fellows. 
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